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Abstract 


A new bivariate function has been found that provides solutions of integrals having 
the form u~’e" du which arise when developing predictions for the behavior of pressure 
within a rigid volume under high vacuum conditions in the presence of venting as well as 
sources characterized by power law transient decay over the range [0,1) for 7 and for u 
> 0. A few properties of the new function are explored in this work. For instance the 7 = 
‘2 case reproduces the Dawson function. In addition, a slight variation of the solution 
technique reproduces the exponential integral for 7 = 1. The technique used to generate 
these functions leads to an approach for solving a more general class of nonlinear 
ordinary differential equations, with the potential for identifying other new functions that 
solve other integrals. 


Introduction 


The evolution of pressure versus time p(t) within a rigid volume V with venting in free 
molecule flow becomes governed by the behavior of gas sources Q contained within the 
volume. Based upon a mass conservation statement for an ideal gas at room 
temperature T [1] 


dp 
V—=Q-C\(p-p,, ). 1 
- (p-p..) (1) 
Venting is described by conductance C to exhaust pressure p,,. For molecular flow 
conditions, C is not a function of pressure, and for this discussion it is assumed p,,= 0. 
Dividing by C and defining characteristic time z = V/Cand reduced gas load q = Q/C, 


7h. 


= q(t)-p. (2) 


Let s=t/r. For sources associated with diffusion-limited outgassing, classical 


theoretical behavior yields qe Js, [2] while for those associated with surface 
desorption qai/s. [3] While these sorts of behaviors are often observed 
experimentally, it appears more generally one must deal with [4] 


q(s)= =, (3) 


where c and 7 are constants for a given system and 0.5 < 7 < 1.0, although this 
parameter range does not represent hard limits. Rearranging and introducing integration 


factor e® 


Ss 


a let) = 0S. (4) 


The solution to this ordinary differential equation for pressure p(t) with some initial 
value p,at t, becomes 


plt)=pie +c 


The matter now becomes one of evaluating the parenthetical expressions. 


Analytical Development 


It appears no general solutions for integrals of the type highlighted in Eq. (5) have 
ever been published, as they are conspicuously absent from references on mathematical 
functions. When 7 is a positive integer their behavior is governed by the exponential 
integral Ei(s), [5-7] and the Dawson function Divs) describes n = *%. [7-9] For this 


integral there also appears to be a relationship with the lower incomplete gamma 
function under certain conditions. [10-12] 


The procedure for direct numerical solution of this integral seems to be quite 
problematic. The integrand has a singularity at s = 0, and for large arguments both the 


numerator and denominator increase without bounds. On the other hand the scaling 
factor ahead of the integral counteracts some of this behavior overall. 


When 0 < 7 <1, let x = %/u, where b =1/(1-77). Then 
S al x 
ea du = be-*’ pe® dx = bG(n, x). (6) 
ou” 0 


It happens that G(v7, x) satisfies its own ordinary differential equation: 


OG head: (7) 
IX 


Properties of Bivariate Function G(7, x) 


The solution of Eq. (6) may be accomplished much more easily numerically using 
Eq. (7) than it can by attacking the integral directly. Eq. (7) may be explored to discover 
a variety of properties of G. 


Integration of Eq. (7) for G across various values of 7 and x produces curves like 
those presented below in Figure 1. 
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FIGURE 1. Function G(, x) for selected values of 7. 


When (n, b, x) = (0, 1, s), G(0, x) reproduces expected limiting behavior of 
s 
e *se°ds ~2(e° - 1)= 1-e%, (8) 
0 


and when 7 = 0.5 the function is equivalent to the Dawson function D(x), with its 
maximum value of about 0.54 at x = 0.92. [7-9] 


As 71, b>, and the function steepens to produce a sawtooth shape with 


infinite slope at x = 1, which appears to anticipate the discontinuity associated with the 
exponential integral. [2] 


Also apparent in Figure 1 is that G’(7, 0) = 1 at x = 0, independent of 7. This 
property is also observed by setting x = 0 in Eq. (7). In terms of physical units, this 
property indicates that regardless of 7» and the fact that the diffusion-desorption 


mechanism starts with a singularity at t = 0, pressure within the volume is initially well- 
behaved with no such influence. 


Less obvious in Figure 1 is the fact that for large values of x independent of 7: 


1 
bxe-1 


G(n, x ~~ 00) > (9) 


In terms of satisfying Eq. (6), recall the solution for the scaled integral was defined by 
bG with s = x”. Substitution of Eq. (9) into Eq. (5) for large s = t/t leads to 


p(t > 2)» — =q(s). (10) 
Ss 


Eq. (10) is the solution to Eq. (2) when p= 0. 
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FIGURE 2. Contour map of function G(7, x). 


A contour map for G(7, x) is presented in Figure 2 essentially spanning the range for 
7 and computed for x out to three in increments of (Av, Ax) = (0.01, 0.01). The linear 
slope in the neighborhood of x = 0 appears as a series of equally-spaced vertical bars, 
and the discontinuity associated with the exponential integral is present as a node at (n, 
x) = (1, 1). The prominent freeform orange cross around x = 0.8 is characteristic of a 
saddle point. Owing to the simplistic integration associated with this figure, one may say 
the saddle point occurs approximately at (7, x) = (0.56, 0.85) with G = 0.54. 


A three-dimensional representation of Figure 2 is depicted below in Figure 3 with the 
n scale reversed to better descry the approach to the singularity at (7, x) = (1, 1). 
Figures 1 & 3 also reveal that another singularity is approached as (n, x) — (0, ~), 
where G(0, « ) > 1, but G(0+,~ ) > 0. 
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FIGURE 3. Three-dimensional contour map of function G(7, x). (7-axis reversed from 
Fig. 2.) 


Scaled Exponential Integral Solution F(y) Properties 


Clearly the change of variable introduced to produce Eq. (6) will not work for 7 = 1, a 
case corresponding to surface desorption. [3] Instead let y =Ins. Then 


S Ql y 
eSs—du=e se” dy =Fly). a) 
go U 
Eq. (11) would generally be described by an exponential integral Ei(&), specifically 
y 
Fly)=e fe® dy = -e Bil. eY ) (12) 


In fashion similar to Eq. (7) it is determined that 


OF eveSd, (13) 
dy 


From Eq. (13) it can be seen that Fily > —co ) > 1, which corresponds to sO. It 
resembles behavior observed for G(7, 0). Also, for large y, analogous to Eq. (9), 


1 
Fly > 0)>eY =—, (14) 
and in physical units 


plt +») >= =d(s). (15) 


Function F(y) is depicted in Figure 4. Its maximum value is approximately 9.60 for 
y ~-2.27,(S ~0.103). 
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FIGURE 4. F(y) behavior featuring function maximum and long-term trending. 


General Solution Approach 


The approach used to identify functions such as F(y) and G(7, x) may be formulated 
in a more general sense beginning with Eq. (2). Casting Eq. (1) in terms of s, the forcing 
term may be introduced as 


q(s) = As)’ (16) 


leading to 


mane (17) 


plt) = pe" 


-s° @" (se her 
[> Fo [> Tay | 


If it is possible to introduce the variable substitution d& = du/f(u) such that u = g(é) 
is a function in terms of the transformed variable, then the scaled integral will be 


du « e92) rege = H(E). (18) 


oO 

b 
on on 
> 
| 
S 
SS 


oe 54, (19) 


If well-behaved, it appears these functions should all rise from some initial 
transformed variable zero with a slope of unity, and for large arguments will approach 
limiting behavior of 


H(& > 2) > 


20 
g'(é) a 


H(&) will also form the core of the solution to integrals of the type arrived at in Eq. (17). 


Concluding Remarks 


A new mathematical function G(7, x) has been discovered that solves integrals 
containing the form u~”e" du through a simple transformation from u to x. The integral 
is solved through use of a property of G as a solution to a particular differential equation. 


This function is useful for describing transient behavior for pressure in vacuum systems, 
and several properties of this function have been explored. 


It has been noted that Dawson’s function is related to certain solutions of the 
diffusion equation used in conduction heat transfer and mass transfer studies, as well as 
certain investigations associated with propagation of electromagnetic radiation. Since 
G(7, X) simplifies to Dawson’s function for 7 = 0.5, it seems logical that the new function 
will offer at least the same level of utility. 


The special case 7 = 1 cannot be solved with the transformation associated with 
development of G, but may be handled by a similar approach. The resulting function 
F(In u) is related to the exponential integral. 


Finally, the approach taken to generate these functions has been generalized to 
identify properties of solutions to integrals containing the form le" /F(u))au. Provided 


certain conditions are met, these solutions will share certain features of short-period and 
long-period behavior. 
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